Eigenvalue statistics of the real Ginibre ensemble 
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The real Ginibre ensemble consists of random NxN matrices formed from i.i.d. standard Gaussian 
entries. By using the method of skew orthogonal polynomials, the general n-point correlations for 
the real eigenvalues, and for the complex eigenvalues, are given as n x n Pfaffians with explicit 
entries. A computationally tractable formula for the cumulative probability density of the largest 
real eigenvalue is presented. This is relevant to May's stability analysis of biological webs. 
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Dyson's three fold way [l[ is a viewpoint on the foun- 
dations of random matrix theory, showing how consid- 
eration of time reversal symmetry leads to three classes 
of ensembles of relevance to quantum mechanics. The 
three ensembles are catalogued by the classes of unitary 
matrices which leave the ensemble invariant — orthog- 
onal (time reversal symmetry is an involution), unitary 
(no time reversal symmetry), and symplectic (time re- 
versal symmetry is an anti- involution) . For an ensemble 
theory of Hermitian matrices, an equivalent characteri- 
zation is that the matrix elements be real, complex and 
real quaternion respectively. 

Both as a concept, and as a calculational tool, the 
three fold way has been highly successful. As a concept, 
allowing for global symmetries in addition to that of time 
reversal gives a classification of the former in terms of the 
ten infinite families of matrix Lie algebras f3| ■ This classi- 
fication now provides theoretical underpinning to funda- 
mental phenomena in mesoscopic physics Q, disordered 
systems [4] , and low energy QCD [5| , in additional to the 
study of the statistical properties of quantum spectra for 
which it was originally intended. A good deal of the suc- 
cess relates to the matrix ensembles of the three fold way 
and its generalization being exactly solvable — analytic 
forms are available for all key statistical quantities, al- 
lowing for quantitative theoretical predictions. 

Soon after the formulation of the three fold way, Gini- 
bre 0] presented as a mathematical extension an analo- 
gous theory of non-Hcrmitian random matrices. The en- 
tries are taken to be either real, complex or real quater- 
nion. Like their Hermitian counterparts, it transpires 
that such random matrices have physical relevance. 

Consider the complex case first. Then the joint eigen- 
value probability density function (PDF) is proportional 
to 
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This can be recognised as the Boltzmann factor for the 
two-dimensional one-component plasma in a disk, or the 
absolute value squared of the wave function for free 
fcrmions in a plane, subject to a perpendicular magnetic 



field and confined to the lowest Landau level Q. In the 
study of chaotic dissipative quantum systems, the statis- 
tical properties of eigenvalues for certain model maps are 
well described by the corresponding statistical properties 
implied by this PDF ||. 

In the case of real quaternion elements, the eigenvalues 
come in complex conjugate pairs. The eigenvalue PDF 
of the eigenvalues in the upper half plane is proportional 
to 
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Up to an extra one body factor YiiLi \ z j ~ %l; the eigen- 
value PDF of the eigenvalues in the upper half plane 
is proportional to the Boltzmann factor for the two- 
dimensional one-component plasma confined to a semi- 
disk, bounded by a dielectric material of dielectric con- 
stant e = along the straight edge 0, [l(| ■ 

Both joint PDFs for the complex and real quaternion 
cases are contained in Ginibre's paper However, in 
the case of real elements, it wasn't until a further twenty- 
five or so years later that the joint distribution was com- 
puted, first by Lehmann and Sommers [XXI ] , then by Edel- 
man [12j. Part of the difficulty is that the joint PDF is 
not absolutely continuous. Rather, there is a non-zero 
probability that for N even (odd) there will be an even 
(odd) number of real eigenvalues for all even (odd) posi- 
tive integers up to N. The final result is that for k real 
eigenvalues (k of the same parity as N), the joint PDF is 
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where A({z p } p= i,..., m ) := n™c/( z * _ z j)- 
(-00,00) while (xj,yj) G R 2 + , M. 2 + := {(x,y) G JR^ : y > 
0}. Integrating over {Xi} Li {xj + iy -} gives the proba- 
bility that there are precisely k eigenvalues. The simplest 
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case in this regard is when k — N (i.e. all eigenvalues 
real), and it is found that the sought probability is equal 
to 2- JV ( A, ~ 1 )/ 4 [12]. For k = 2 an evaluation in terms of 
a single definite integral has been given recently in [l3j |. 
while [HI reduces the calculation for general k down to 
an expression of the same computational complexity as 
our Eq. (fl"Tj) below. 

Perhaps the first applied study to draw attention to 
the eigenvalues of random real matrices was that of 
May [l5j], in the context of the stability of large biolog- 
ical webs. The very general setting [16|] is to consider 
an A-dimensional vector x(t) with components speci- 
fied as the solution of the coupled first order system 
dxi(t)/dt = Fi(x(t)), (i = 1,...,N) for some nonlinear 
functions i^. Assuming an isolated equilibrium solution 
x° , linearization about this point leads to the linear ma- 



trix differential equation 

dm 

dt 



Ay(t) 



(4) 



where A is an N x N matrix. The system is stable if 
all eigenvalues of A have a negative real part. To model 
the effect of random coupling between components on a 
stable equilibrium, the matrix A is written A = —In + B 
where B is a dilute matrix (fraction 1 — c of its elements 
zero) with mean zero and variance s 2 . The May crite- 
rion asserts that stability requires sV Nc < 1. Indeed in 
the case c = 1 this is consistent with limit theorems for 
the spectral radius of random real matrices proved sub- 



sequently [17j, ll8| . Neural networks are further examples 
of complex webs to which such a random matrix based 
stability analysis is relevant 1|| 2(| 21 1 . The results ob- 
tained below allow the evaluation of the probability of 
stability in the borderline case of the May stability cri- 
terion, 



'Nc = 1. 



(5) 



As with the matrix ensembles of Dyson's three fold 
way, all correlations and a number of key distributions 
for the complex and real quaternion Ginibre ensembles 
have been calculated exactly 0, [23| ■ The Fourier trans- 
form of the two-point correlation (structure function) is 
a quantity of key importance to the plasma and fermion 
interpretation of |T]), while the decay of the two-point 
function along the boundaries indicates general physi- 
cal principles (non-zero dipole moment of the screening 
cloud in the case of fl}; vanishing dipole moment for 
©). Further, the distribut ion function for the spacing 
between eigenvalues in the bulk can be compared against 
data obtained from dissipative maps [H, while the den- 
sity fluctuations in a large disk within the bulk indicate 
further general physical principles [221 ] . 

In contrast to the situation for (fl]) and ([2|), the corre- 
lations and distributions have not in general been com- 
puted for the real Ginibre ensemble. Exceptions are the 



density of real eigenvalues [24| 



2tt\ T(N-l) 

with T(p,x) := f™t p - 1 e- t dt, ~/{p,x) :=J* i^e"* dt, 
and the density of complex eigenvalues [l2| 

^2T(N-l,X 2 + y 2 ) 2 , 
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Further, with Z k ^ N _ k y 2 [u, u ] denoting the canonical av- 
erage of nf=i u i^i) llj=r fc ^ 2 u ( x j + w ^h respect to 
Q , it has been shown in [25[ (taking N even for definite- 
ness) that 

JY 

Z N [u, u] := y]Z ki ( N _ k) /2[u,u\ 

0-JV(JV+l)/4 

= ttJV rfJ/ J >i l a oA U ) + Pj,k( U )kk=l,...,N, (8) 

where, with pi(x) an arbitrary monic degree I polynomial 
and z := x + iy, 



&j,k(u) = / dxu{x) I dyu(y) 
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Pj,k(u)= 2i dxdyu(z)e y x erfc(v / 2y) 

x(Pj-i(z)p k -i(z) - p k -i(z)pj-i(z)). (10) 

It is the purpose of this Letter to report that all the re- 
sults ©-([SI) can be generalized, thereby fully exhibiting 
the real Ginibre ensemble as exactly solvable. For con- 
venience it will be assumed throughout that N is even. 
We first observe that with the second argument u on the 
LHS of © replaced by an arbitrary function v = v(x,y) 
the equality remains valid with u in /3j tk (u) replaced by 
v. With PN2n denoting the probability that 2n out of 
the N eigenvalues are real, it then follows by choosing 
v = 1 , u = £ that 
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(cf. Eq. (11) of Ref. [l4). 

As is well known in random matrix theory [ij [23| , the 
correlations of a Pfaffian generating functional (JSj> are 
themselves Pfaffians. However in general this form in- 
volves the inverse of the matrix in (151) with u = 1. To 
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make this explicit, one seeks to choose the polynomials 
{pi(x)} to have the skew orthogonality property 

Qt2j,2fc(l)+/32i,2fc(l)=Q!2j-l,2fc-l(l) + /92j-l,2fc-l(l) = 0, 
Q!2i-l,2fc(l) + /32i-l,2fc(l) = Tj-i8j t k- (12) 

Our key result is that the very simple choice 

p 2 , (z) = x 2 \ P2 j+ i(x) = - 2jx 2j -\ 

rj-x = 2v / 2^r(2j - 1) (13) 

validates (fT2|) . With this established, and q2j(z) := 
— p2j+i(z), q2j+i(z) :— p2j(z), one finds for the corre- 
lations between complex eigenvalues 
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where S c (w,z) := J2j=i Pj-i{ w )<lj-i ( z )/ r [(j-i)/2] and 
zj := Xj + iyj. In the case n = 1, the Pfaffian equals 
5 c (zi,zi) and ([7]) is reclaimed. In the case n — 2 the 
Pfaffian equals S^zi, zi)S c (z 2 , z 2 )+S°(zi, z 2 )S c (zi, z 2 )— 
S c (z 1 ,z 2 )S c (z 1 ,z 2 ). 

Similarly, the correlations between real eigenvalues are 
computed as 



p\ n) {x u ...,x n ) 
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with S*(x,y) = ^r(x,y), D*(x,y) = ^S T (x,y) and 
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In the case n = 1 this gives p/nix) = S T (x,x) and (JB|) 
is reclaimed. In the limit N — ► 00 with x,y fixed (p~6|) 
simplifies to 



/ r (^,2/) 



(17) 



implying the correlations decay at a Gaussian rate. In- 
tegrating (O over x S (—00,00) gives the mean number 
En of real eigenvalues, which is computed [24[ to have 
the large N asymptotic form v / 2N/tt(1 - 3/87V -...). 
The variance Vn of this same number is computed in 
terms of the two-point correlation according to Vjv = 
I-oo dx J r Zo d V : Pp)( x >v) + e n, Pp)(x,y) ■= p\ 2) (x,y) - 
P^{x)p r ^(y). We read off from (|T5j) that 

p T ^)(xi,x 2 ) = -S r (xi,x 2 )S r (x 2 ,xi) 

-(sgn(xi - x 2 ) + T (xi, x 2 )j D 1 (x!, x 2 ). (18) 



The quantity p^.(x,y)/p^(x) is integrable in y showing 
that for large N, 

V N ~E N (l + limjl/p {1) (0)) J p\l ) {Q,y)dy) 
= (2-V2)E N , (19) 

with the final equality making use of ([IB"]) and (fiTf . 

We draw attention to quantitatively similar results 
which hold for the zeros of the random polynomial p(z) = 
1/2 



v-^iv ( n\ 
^3=0 \ j ) 



otjZ J , where the 



.d. real Gaussian 



random variables. This has the interpretation in quan- 
tum mechanics as a random superposition of states with 
spin N/2. Moreover, the function p(e 1 ^ cot 9/2) vanishes 



at the values of 



) on the sphere corresponding to 



the stereographic projection of the zeros Zj of p(z), giv- 
ing the Majorana parametrization [28| . The analogy with 
the present problem is that the mean number of real ze- 
ros is proportional to y/N, as is the variance, and the 
correlations decay as Gaussians [29[ . A distinction is the 
lack of a boundary for the eigenvalue support, which is 
distributed as a Cauchy distribution. 

We remark too that although not reported on here, 
the correlations between real and complex eigenvalues 
can be written as a Pfaffian. Furthermore, we antici- 
pate that the par tially symmetric real Ginibre ensemble, 
introduced in [llj, will also yield to the present strategy. 

To leading order the support of the eigenvalue densi- 
ties ©, (J7J| is the disk \z\ = \/~N, as is consistent with 
the formula for the boundary of the May stability cri- 
terion (here c = 1 and s = 1; however the variable s can 
be reinstated by scaling z 1— » z/s throughout). Setting 
x = \/N + X and taking N — > 00 in ([6]) gives for the 
limiting edge profile 



1 /l 



(-(l-erfV2X) 
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(1 + erfX) 



(20) 

For any fixed angle away from the real axis, as N — > 00 
the density of complex eigenvalues near the boundary of 
support is radially symmetric, and the same as that in 
the complex Ginibre ensemble. Writing the radius r as 
r = \fN + R, for N — > 00 this has the explicit form 



pl ) (R) = -^{l-eiiV2R), 



(21) 



and is thus equal to twice the first term in (p0| . 

Suppose now that the variance of the Gaussian en- 
tries is reinstated as the variable s 2 . Let R r be the 
event there are no real eigenvalues, or all real eigen- 
values are less than r. By scaling of the eigenvalues, 
Pr(R s ^ N - +sr ) is independent of s, and for N — > 00 it 
is an order 1 function of r. The latter can be writ- 
ten as an infinite sum over the limiting n-point edge 
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correlations, or equivalently as a Fredholm determinant 
of the integral operator with kernel given by the edge 
limit of the general entry in (|15[) . For r large one 
has limAr^oo Pr(i? sv ^ +sj ,) ~ 1 - J*™ (X) dX, show- 
ing that the corresponding PDF decays as a Gaussian. 
For general N a practical formula for computing this 
probability is in terms of the generating functional l|8|). 

PT ( R aVN+sr) = Z N[Xxe(-oo,sVN+sry 1 ]' whcrC XA = 1 

if A is true, \A — otherwise, and with the polyno- 
mials in |9j), (fTTj]) chosen according to (Tl3|) . Numeri- 
cal values of Pr(ii g so computed are tabulated in 
Table |T] for successive even values of N. The quantity 

Pr (-^ s ^) := ( Pr (- R sv / lv) _ Pat,o)/(1 - Pjv.o), also listed 
in Table U gives the probability that all real eigenvalues 
are less than syN, given that there is at least one real 
eigenvalue. In the case s 2 — l/N, this corresponds to 
the probability that all non-oscillatory solutions of the 
linear system §4§ are stable, given that there is at least 
one non-oscillatory solution. 



N 






2 


0.81444 


0.737579 


4 


0.793864 


0.756706 


6 


0.784485 


0.762255 


8 


0.778838 


0.764193 


10 


0.774963 


0.76475 


12 


0.772092 


0.76469 


14 


0.769855 


0.76434 


16 


0.768048 


0.76385 



TABLE I: Tabulation of two probabilities, specified in the 
text, relating to the probability that all real eigenvalues of an 
N x N Gaussian real matrix, entries of mean zero, variance 
s 2 , are less that sy/N. 



With this study, building on the contributions of 
Lehmann and Sommers [ill ]. Edelman [12J, Kanzieper 
and Akemann [3] and Sinclair (2||, the problem began 
by Ginibre over forty years ago of calculating the statis- 
tical properties of the eigenvalues of non-Hermitian real 
Gaussian matrices is solved. As a consequence the dis- 
tribution of the largest real eigenvalue is presented in 



a computable form. The largest real eigenvalue deter- 
mines the stability of non-oscillatory solutions in May's 
(lit analysis of biological webs. 

The work of P JF has been supported by the Australian 
Research Council. 
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